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ABSTRACT 

Limits to the Gunn-Peterson effect due to neutral hydrogen have generally been obtained by 
modeling the observed Lya forest as a superposition of absorption lines with Voigt profiles arising 
from clouds of photoionized gas, and a hypothesized uniform continuum arising from the intergalactic 
medium. However, owing to the formation of structure by gravitational instability, a photoionized 
intergalactic medium should be inhomogeneous on scales larger than the Jeans scale, and therefore 
the optical depth should fluctuate. Such a fluctuating continuum can always be modeled as a 
superposition of lines, but this decomposition does not necessarily have a direct physical meaning. 

We present a calculation of the evolution of the density in voids in a photoionized intergalactic 
medium, using the Zel’dovich approximation and another analytical approximation which we argue 
should be more accurate in this regime. From this, we calculate the probability distribution of 
the Gunn-Peterson optical depth in terms of the amplitude of the primordial density fluctuations. 
Over most wavelengths in a quasar spectrum, the optical depth originates from gas in underdense 
regions, or voids. Individual absorption lines should be associated with overdense regions, which we 
do not treat here. This causes the median Gunn-Peterson absorption to be lower than the value 
for a uniform medium containing all the baryons in the universe by a large factor, which increases 
as gravitational collapse proceeds. The Gunn-Peterson effect is the only known method to directly 
observe underdense matter in the universe, and it can be sensitive to the primordial fluctuations 
even in the non-linear regime. In particular, in the He II Gunn-Peterson effect recently detected by 
Jakobsen et al. , gaps in the absorption are a very sensitive probe to the most underdense voids. 

We apply our calculations to the observations of the intensity distribution in a 2 = 4.11 quasar by 
Webb and coworkers. We show that if Lya clouds arise from gravitational collapse, their observations 
must be interpreted as the first detection of the fluctuating Gunn-Peterson effect, with a median value 
Top — 0.06 at z = 4. If the linearly extrapolated rms density fluctuation at the Jeans scale for the 
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photoionized gas were close to unity at this redshift (which is the case in typical low-density models 
with cold dark matter), then tgp should be ~ 1/5 of the optical depth that would be produced by a 
uniform intergalactic medium. This is consistent with the predicted baryon density from primordial 
nucleosynthesis, and the intensity of the ionizing background derived from the proximity effect. 
From the numerical simulations of Cen et al. , such models also predict correctly the number of 
Lya absorption lines observed. For theories with much larger density fluctuations (such as standard 
cold dark matter), we argue that, given the observed number of lines with Nhi ^ 10^'*cm“^, the 
Gunn-Peterson optical depth should be much lower than observed; this needs to be investigated in 
more detail using numerical simulations. 

Subject headings: intergalactic medium - quasars: absorption lines - cosmology: large-scale structure 
of Universe 
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1. INTRODUCTION 


If the baryonic material that is seen today in galaxies (and X-ray gas) had in the past been in 
the form of neutral gas distributed in the intergalactic medium, a trough should have been seen in 
the spectrum of any sufficiently distant quasar lying behind such a neutral medium at wavelengths 
below the Lya line. The trough is caused by scattering of photons as they are redshifted through 
the frequency of the Lya line (Gunn & Peterson 1965). The absence of the trough implies that 
the medium had been ionized before the epoch of the highest redshift quasar known [z = 4.89; 
Schneider, Schmidt, & Gunn 1991). 

If the ionization is due to photons, an abundance of neutral gas in photoionization equilibrium 
should still be present. This should yield a depression of the flux below the Lya line, which would 
be constant for homogeneously distributed neutral hydrogen. What is observed instead is the Lya 
forest, caused by a rapidly varying optical depth to Lya scattering with wavelength. This reflects a 
highly inhomogeneous distribution of neutral hydrogen. 

Since the discovery of the Lya forest, there have been several attempts to detect a uniform 
Gunn-Peterson effect, by subtracting the contribution to the average depression of the flux of the 
individually detected absorption lines. This has led to progressively more stringent lower limits to 
the intensity of the ionizing background, if the density of the intergalactic medium is assumed to 
be close to the baryonic density required by primordial nucleosynthesis (Steidel & Sargent 1987; 
Jenkins & Ostriker 1991; Webb et al. 1992; Giallongo et al. 1993, 1994). Since this intensity is also 
constrained by the observation of the proximity effect (e.g., Bajtlik, Duncan, & Ostriker 1988), this 
could eventually lead to the conclusion that the intergalactic medium must be collisionally ionized. 

These limits are based on the assumption of the presence of a uniform intercloud medium, as 
distinct from individual clouds producing the absorption lines. However, there does not seem to 
be any theoretical justification for such a distinction for the case of a photoionized intergalactic 
medium. The galaxies and clusters we observe today tell us that the universe is not homogeneous. 
They must have collapsed from density fluctuations which were already substantial even at the 
highest redshifts of observed quasars. Moreover, in any theory of hierarchical gravitational collapse, 
fluctuations on subgalactic scales are even larger, and collapse at an earlier epoch. Since the Jeans 
mass for a photoionized medium (with a temperature T ~ 10'* K) is ~ 10*°[(1 + 2:)/5]“^/^ Mq, the 
gas down to these very small scales should have a fluctuating density, and a uniform intergalactic 
medium cannot exist. Therefore, the Gunn-Peterson effect cannot be uniform. 

Primordial fluctuations on small scales should cause the photoionized gas to collapse first in 
objects of the Jeans mass at T ~ 10'* K. This idea led to the minihalo model (Rees 1986; Ikeuchi 
1986), and other alternatives also based on gravitational collapse as the origin of the Lya clouds 
responsible for the observed absorption lines (e.g.. Bond, Szalay, & Silk 1988). Any residual Lya 
absorption between the lines would then mostly come from underdense regions, or voids, that are 
left as the gas collapses gravitationally and forms structures. The optical depth in such spectral 
regions can be reduced by a large amount if voids are underdense by large factors. Yet, because the 
pressure of the photoionized gas prevents it from collapsing on scales smaller than the Jeans scale, 
an intergalactic medium pervading all space must be present, and the inhomogeneous distribution of 
gas in voids must yield a fluctuating Gunn-Peterson effect. 
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Recently, Cen et al. (1994) presented a cosmological simulation of the Lya forest, computing the 
evolution of photoionized gas when structure on the relevant scales of 50 — 500/i“^kpc in comoving 
units collapses gravitationally. They showed how the high contrast of the lines in the Lya forest can 
arise naturally in structures resembling Zel’dovich pancakes (Zel’dovich 1970; Sunyaev & Zel’dovich 
1972). The contrast arises from the overdensity of the gas, and the fact that the neutral density is 
proportional to the square of the total density. 

In this paper, we use the Zel’dovich approximation and an ad hoc variant of it to calculate the 
probability distribution function of the Gunn-Peterson optical depth. This semi-analytical method 
can be useful for underdense regions where the gas has not yet been shocked in any collapse. In § 2 
we describe our method, and in § 3 we present the results for various models of structure formation. 
§ 4 shows how our results can be used to interpret the observations of the weakest lines in the Lya 
forest as a fluctuating Gunn-Peterson effect, and discuss future applications of these observations 
as a probe to the primordial density fluctuations on small scales. Finally, in § 5 we summarize our 
results. 


2. METHOD 


In this section, we present the method by which we estimate the distribution of the optical depth 
to Lya scattering caused by different regions in an inhomogeneous intergalactic medium, where the 
gas is falling into gravitational potential wells arising from initially gaussian fluctuations in the 
cold dark matter (GDM) distribution that grow by gravitational instability. We use the Zel’dovich 
approximation, and also what we call the “modified Zel’dovich approximation”, by which we attempt 
to improve the prediction for the evolution of the density in underdense regions, or voids, as the rms 
fluctuation reaches values of order unity. In §2.1, the relation between the optical depth observed 
in a quasar spectrum and the density and velocity fields in the intergalactic medium is presented. 
§ 2.2 shows the relation between the GDM and gas density fluctuations as affected by the finite gas 
pressure. In §2.3, we discuss the usual Zel’dovich approximation and our modification; and in §2.4, 
we describe the formalism to calculate the distribution of the density and the optical depth with 
these approximations, for gaussian initial density fluctuations. 

2.1. The fluctuating Gunn-Peterson optical depth 
of an inhomogeneous intergalactic medium 

Throughout this paper, we assume that the gas in the intergalactic medium is photoionized and 
in ionization equilibrium with an ionizing radiation background of uniform intensity. Let the gas 
density at any point in space, with comoving coordinates x, be Pg{'x.) = pg Vg{pf), and let the density 
of neutral hydrogen be Pn{^) = Pnu i^n(x), where pg is the average gas density in the universe, and 
Pnu is the density of neutral hydrogen for a uniform gas of density pg in ionization equilibrium 
with an ionizing background of the same intensity as in the nonuniform case. (For example, in an 
isothermal and mostly ionized medium, = I'g-) The neutral hydrogen yields an optical depth 
due to Lya scattering along the line-of-sight to a quasar which can be observed as a function of 
wavelength. The wavelength is related to the component of the total velocity of the gas along the 
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line of sight, y, by A = Anya (1 + -Zo) (1 + ^jc), where Anya is the Lya wavelength, and the velocity 
obeys <C c, and is measured with respect to an arbitrary origin at redshift zq. 

For a uniform medium, this optical depth is given by (Gunn & Peterson 1965) 

Tu = 4.65 X 10® h{pnu/Pg) (1 + zT[Ho/II{z)] , (1) 

where Ug is the gas density divided by the critical density, Hq = 100/ikm s“^ Mpc~^ is the Hubble 
constant at present, and H[z) is the Hubble constant at redshift 2 . In the general case of a medium 
with varying density, the optical depth contributed by any point is increased proportionally to the 
neutral hydrogen density, and inversely proportionally to the gradient of the total velocity along the 
line of sight of observation. In general, there can be more than one point along the line of sight 
having a velocity V. Let these points have a comoving coordinate where x is the component 

along the line of sight. Ignoring the effects of thermal broadening, we then have 


t{V) = Tu 




a 

\dVjdx 


(2) 


where a is the scale factor and dots denote derivatives with respect to cosmic time t, so that 
H[z) = d/a. When the density fluctuations are small, there is a unique position to every velocity, 
but as the fluctuations grow velocity caustics will form, and within those x[V) is multivalued. 
Velocity caustics generally appear earlier than caustics in real space (e.g., McGill 1990). 

The derivative of the velocity appearing in equation (2) can be written as dV/dx = 
[dVi/dxj) ei Cj, where is the unit vector along the line of sight. The total velocity is V = dx + v, 
where v is the peculiar velocity. Thus, 

dV (. dvi \ 


2.2. Growth of density perturbations in a gas with finite pressure 

In order to encompass a number of the currently popular models for the growth of structure 
in the Universe, we allow for cold (GDM) and hot (HDM) dark matter components. The GDM 
has a negligible temperature and can grow perturbations on all scales of interest, which produce 
gravitational potential wells that affect the dynamics of the gas. The HDM, on the other hand, is 
so hot that its Jeans scale is much larger than the scales of interest, and it will be taken as uniform. 
In the present section, we can ignore the presence of the HDM. 

The temperature of the gas in the intergalactic medium is sensitive to the history of reionization, 
since at low redshifts the thermal equilibrium time is longer than the Hubble time (e.g., Rees & Setti 
1970). Assuming that the dominant heating and cooling mechanisms of the intergalactic gas are, 
respectively, photoionization heating and adiabatic cooling, Miralda-Escude & Rees (1994) found 
that, for redshifts 2^5 and for the average baryon density, the gas temperature will tend to a 
value T sa To(l + z)^ within a Hubble time, where To sa 3000 K and depends on the spectrum of the 
radiation, and /3 sa 0.9. For other densities, the temperature is T oc (see their Figs. 2 and 3). 
For mathematical convenience (see next paragraph), we take (3=1. Since the photoionization rate 
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is oc pn, and the recombination rate is cxi (assuming <C Pg), the equilibrium density of 

neutral gas obeys i^n ~ ^g'^- 

In the linear regime (i.e., for small perturbations to a uniform-density background), and if the 
temperature decreases as T = To(l -f -z), then the power spectrum of the gas density fluctuations is 
related to that of the cold dark matter by 


PM 


PcDM{k) 

[i + {k/kjy]^’ 


(4) 


with the Jeans wavenumber 


1/2 

, _ f ‘iT^tJ'mpGpcDM,o \ 

“ 1, jksTo ) ' ^ ^ 

Here, p is the mean molecular weight in atomic units (= 0.6 for fully ionized hydrogen and helium 
in cosmological abundances), rrip is the proton mass, G is the constant of gravitation, pcDMfi is the 
current average density of the cold dark matter, 7 is the value of the derivative din p/d In p in the 
perturbations (7 = 4/3 under the assumptions made above), and ks is Boltzmann’s constant. A 
sketch of the derivation of this relation for a flat, CDM-dominated universe can be found in Peebles 
(1993, pp. 635-636), and the extension to the present, more general case is relatively straightforward. 
When Pcdm(^) increases less rapidly than k as k ^ 00 (which is expected to be the case because 
of smoothing processes before recombination), the velocity dispersion of the cold dark matter in 
collapsed objects decreases toward small scales; when this velocity dispersion drops below the sound 
speed of the gas, pressure prevents the gas from collapsing in the dark matter potential wells. The 
gas is therefore smooth on scales smaller than the Jeans mass, and the rms fractional fluctuations 
of the gas density have a finite value, which we define as aj. 

In addition to the small-scale smoothing that appears already in linear perturbation theory, 
the pressure also affects the nonlinear evolution of the gas, and this can be important on scales not 
much larger than the Jeans scale. These effects cannot be included in our treatment. 


2.3. Gravitational evolution 


In this section, we study the gravitational evolution of the gas, treated as CDM with the power 
spectrum modified as in equation (4). Our treatment is for a general Friedmann-Robertson-Walker 
universe with arbitrary values of the density parameter fi and the cosmological constant A. We take 
the CDM to have a fraction fcDM of the total mass, with the remainder being a uniform, static 
background (HDM). 

It is useful to introduce the velocity variable u = x' = v/(od), where x are the comoving 
Eulerian coordinates, the primes denote derivatives with respect to the expansion factor o, and v 
is the peculiar velocity. The deformation tensor Vu = 5u/5x can be decomposed in the standard 
way. 


V iUj — ^ OSij -f “t“ LOij , 


(6) 
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where 0 = V-u is the trace, and the last two terms are, respectively, the traceless symmetric and 
antisymmetric parts of the tensor. The evolution of the normalized density \v\ = pjp (the variable 
V diverges and changes sign when a caustic is crossed) the velocity divergence 6, the shear, and the 
vorticity are determined by 


ly' = —6v, 


( 7 ) 


, 2-0 , 

e' = - -e - 

3 a 


- ^fcDM{v — 1 ), 


(8) 




2 g 2 

a o 

3^^ ^ a {dx,dx, 3^ ^ ’ 


V I ! 


and 
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iJij — T^ikUlkj — ^ik^kj- 


Here, S — ^ to — uj'ijuj'ij^ 


p{a) 


fio/2 

fio T (1 — f2o — nA,o)® + 


( 9 ) 


( 10 ) 


( 11 ) 


and 


q{a) 


adj Qo/2 — f^A 0^^ 

fio T (1 — f2o — nA,o)® + 


( 12 ) 


is the cosmological deceleration parameter, given in terms of the present values (denoted, here and 
below, by a subscript 0 ) of the fractions of the closure density contributed by matter, Q, and by a 
cosmological constant, Qa = A/3H^ (e.g., Peebles 1993, p. 100). The gravitational potential p is 
the solution of the Poisson equation 


— fcDMiy — l)/a. (13) 

The only term in equations (7) through (10) which does not only depend on the local values of 
the density and the spatial derivatives of the velocity field is the term proportional to d^p/dxidxj 
in eq. (9). This term makes it impossible to make an exact, local calculation of these quantities in 
terms of their initial values (Kofman & Pogosyan 1994). However, a few local approximations to 
these equations have been proposed and used (e.g., Zel’dovich 1970; Bertschinger & Jain 1994), with 
variable success. In this paper, we use the Zel’dovich approximation (Zel’dovich 1970; Shandarin & 
Zel’dovich 1989) and a variant of it, both of which are described below. All local approximations 
necessarily break down at orbit crossing. 
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The growing modes in linear theory have no vorticity, and if Wij = 0 initially, it will be true 
until orbit crossing occurs (Peebles 1980, §22 and 23). Thus, we set Wij = 0 in our approximations. 

We write the Lagrangian coordinates of a given mass element as y, defined so as to coincide 
with its comoving Eulerian coordinates x at the initial time, o sa 0. Now, in both approximations 
we make the assumption that 


5u , 

— [y,a) = D{y,a) — {y,a 


0) 


(14) 


i.e., that the time-evolution only multiplies the Lagrangian-space deformation tensor by a (possibly 
position-dependent) scalar D'[y, a). The Jacobian matrix of the transformation from Lagrangian to 
Eulerian coordinates is 





OU' 

6ij + D{y, o)^(y, 0) = [1 - D{y, o)Ai] 6ij 

oyj 


(15) 


where the last expression is valid in a local Cartesian coordinate system with axes aligned with the 
eigenvectors of the initial deformation tensor, and defines the eigenvalues A^. The determinant of 
this matrix has to be 


= (l-T)Al)(l-T)A2)(l-T)A3) = ^^-^ (16) 

in order to satisfy the continuity equation (7). The Eulerian deformation tensor becomes 

duj _ duj dyk _ D'\ 

dxj dyk dxj 1 — D\i 

The dynamical equation (eq. [ 8 ]) can be used to determine the evolution of D. Writing all 
variables in terms of D, and defining ? 7 i = Ai -f A 2 -f A 3 , 772 = A 1 A 2 -f A 1 A 3 -f A 2 A 3 , and 773 = A 1 A 2 A 3 , 
we have 


(17) 



D" I ^ ~ g jj/ ^ ^fcDMP VlD - V2D^ + V3D^ 
a 771 - 2 r] 2 D + SrjsD^ 

This equation, when linearized in D, becomes independent of the initial eigenvalues and 
reproduces the familiar Zel’dovich approximation (hereafter ZA; Zel’dovich 1970; Shandarin & 
Zel’dovich 1989). This is known to be the correct linear limit of the full Lagrangian dynamical 
equations (note that the problematic shear term does not contribute to this order), and gives the 
exact dynamics in the one-dimensional case. In this approximation, the full dynamics is described by 
a single growth factor D[a), the same at all points of space. It gives a remarkably good description 
of the gravitational evolution of matter up to (but not beyond) shell crossing. 

In our modified Zel’dovich approximation (MZA), we solve the nonlinear equation (18) in order 
to obtain D[\,a) independently for each set of eigenvalues A = (Ai,A 2 ,A 3 ). This approximation 



gives the exact evolution as long as the velocity field and gravitational potential satisfy the tensor 
relation (obtained from evaluating eq. [9] in terms of D and its derivatives, and comparing to eq. 
[ 18 ]) 


Vu = a(x, o)VV^, (19) 

where a(x, o) is an arbitrary scalar function. This relation is satisfied by the growing modes in 
the linear regime (where both this approximation and the ZA reduce to the exact linear theory), 
but does not remain correct in the general case. However, it is always satisfied in some simple 
geometries, in particular for planar, cylindrical, and spherical symmetry. Thus, in these geometries 
the approximation gives exact results up to shell crossing. This makes it clearly better than the 
ZA, which only gives exact results in a planar symmetry. In particular, in a spherical void in an 
Einstein-de Sitter cosmology (p = g = 1/2) with fcDM = 1; the ZA gives an asymptotic density 
decrease v oc instead of the correct result (reproduced by the MZA), v oc In fact, in 

a general void, if the tidal forces from collapsed objects near the edges of the void are neglected, 
particles should simply keep their original velocities in the limit where the density in the void is 
very low, and therefore the shape of the deformation tensor in Lagrangian coordinates should stay 
constant, which is the approximation (14). This suggests that, while the ZA generally exaggerates 
the emptiness of the voids, the MZA may describe them more accurately. 

However, in addition to the inconvenient fact that the growth factor depends on the initial 
condition, the MZA has the problem that in cases when ? 7 i < 0 but not all three eigenvalues have 
the same sign, the integration will reach a singularity (corresponding to a zero of the denominator of 
the term on the right-hand side of eq. [18]) where D' diverges, and beyond which the solution is not 
single-valued. Physically, the density initially decreases but reaches a minimum at which it turns 
around because the matter is contracting in at least one direction, finally leading to collapse in the 
ZA. In the MZA, the shear diverges at the turn-around point, leading to an unphysical singularity. 
In order to avoid this complication, we limit ourselves to cases where all < 0 (expansion in all 
directions) when using this approximation. 

The MZA is similar to the approximation of Bertschinger & Jain (1994) in that it leaves 
the principal axes of the deformation tensor of a given Lagrangian matter element fixed in space. 
However, it differs in that the Bertschinger-Jain approximation predicts collapse into filamentary 
structures, whereas the MZA clearly produces “Zel’dovich pancakes.” 

Finally, we determine the initial conditions for equation (18). In the limit o —> 0, one has 
D oc o^, with ^ = [(1 -|- 24/c'dm)^^^ ~ l]/4 (Bond, Efstathiou, & Silk 1980). The value of the 
constant of proportionality in this relation is obtained as follows. We take {ri\) = 1, so the rms 
fractional density fluctuation of the gas at small scales (as calculated in linear perturbation theory) 
should be aj = D{a) as obtained in the ZA. Thus, one can use the modified power spectrum from 
equation (4), together with the (directly or indirectly) measured amplitude at some other scale (e.g., 
the popular as, the rms fractional density fluctuation in a sphere of radius 8/i“^ Mpc, again in linear 
theory) to infer the present value oi D = aj in the ZA, fixing the normalization. In the MZA, we 
simply require the asymptotic behavior of D[y, o) as o —> 0 to be the same as that of D[a) in the 
ZA. 


2 . 4 . Probability distribution functions 
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Doroshkevich (1970) calculated the PDF of the initial eigenvalues of the deformation tensor 
for a randomly chosen point in Lagrangian space for random Gaussian fields. We use this PDF to 
generate numerically a random set of values of the eigenvalues, for which we calculate the values of 
p and T expected in different models at different redshifts (as explained in §§2.1 - 2.3), and finally 
combine the obtained values to obtain PDFs for these two variables. 

In order to generate random values for the eigenvalues A^, we find it convenient to change to 
new variables R, S, and Z by 


= 3 + 7! ^ y3(r^]} - 

1 r 2RS] , , 

(Bertschinger & Jain 1994 use a similar change of variables.) In order to have Ai > A 2 > A 3 , we 
require 0 < JZ < 00, —1/2 < S < 1/2, and —00 < Z < 00. Unlike the eigenvalues A^, these new 
variables are independent of each other, with joint probability distribution 


which makes it easy to generate random values by standard algorithms (Press et al. 1992). 

Kofman et al. (1994) point out that, in general, the density in Eulerian space is the result of 
adding the densities contributed by several particles (“streams”) coming from different Lagrangian 
coordinates. This occurs after the first “pancakes” have collapsed. In the scenario of interest to us, 
the gas is shocked when that takes place, invalidating the approximations of §2.3 in these regions. 
The single-stream probability distribution of the eigenvalues in Eulerian space is Q(Ai, A 2 , A 3 )/|i^|. 
Thus, the single-stream density PDF as sampled in Eulerian space takes the form 


v{p) = J dAi dX2 dX 3 (p _ -p\^\) (22) 

(Kofman et al. 1994). 

In order to obtain a PDF of the optical depth that can be compared with observations of quasar 
spectra, one must take into account that the latter do not sample the gas uniformly in Eulerian space, 
but rather in velocity, along the line of sight. It is convenient to rewrite eq. (3) as dV/dx = dp, 
where (in Cartesian coordinates aligned with the principal axes of the gas element of interest) 


M = 


dij 


dui\ d\nDDXisl 

— eiSi = 1 --- > -— 

dxj J d In o 1 — DXi 


The optical depth contributed by one mass element is 


(23) 


ImI 


r = 


(24) 
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Now, we calculate the PDF of the optical depth in an “observational space” with velocity as the 
coordinate along the line of sight, and Eulerian distance in the transverse directions. In this space, 
the single-stream PDF of the eigenvalues for a fixed orientation of the line of sight with respect to 
the eigenvectors of the deformation tensor is |/x/i^|Q(Ai, A 2 , A 3 ), and the single-stream PDF of the 
optical depth can be written as 


/ 


T{t) = / d\id\2d\3 


sin 6d6 d<f) 
47r 




iMl 


(25) 


Here, 6 and (f) are angles giving the direction of the unit vector e^. The single-stream PDFs are equal 
to the total (observable, multi-stream) PDFs, only if the contribution from regions which have gone 
through caustics in the respective space is negligible. In this case, it can be shown that the integral 
Ns of the PDF over all possible parameter values is unity. In general, NTj > 1 is the average number 
of regions contributing to the optical depth at a given point in the spectrum. 

As an estimate of the low-density end of V{p) and the low-optical depth end of T{t) in the 
general case, we use eqs. (22) and (25), except that we integrate only over those values of the 
integration variables which correspond to regions that have not yet gone through caustics in the 
space of interest, i.e., for which 1 — DAi > 0 in the case of V[p), and for which both 1 — DAi > 0 
and /X > 0 for T{t). For small p and r, this should be a better approximation than integrating 
over the whole parameter space, since in the latter case we would also include regions which in our 
formalism have gone through a caustic and are reexpanding, but in reality are shocked as soon as 
they approach the caustic, and do not expand afterwards. Of course, the high-density regions are 
only shocked once they reach caustics in real (Eulerian) space rather than the “observational space” 
used to calculate T{t), but this should not make a difference for the behavior at small r. 


3. RESULTS 


We use the procedure described in §2 to obtain PDFs for four different models of structure 
formation: 1) standard CDM (sCDM) with Qq = 1 and h = 0.5; 2) mixed dark matter (MDM) with 
Clo = 1, h = 0.5, and fcDM = 0.7; 3) CDM in an open universe (CDMo) with Qq = 0.3, Qao = 0, 
and h = 0.65; and 4) CDM in a low density, flat universe (CDMA) with Qq = 0.4, Qao = 0.6, and 
h = 0.65. 

In order to normalize the growth factor in each of these models, we take the present value of 
the linearly extrapolated rms fractional mass fluctuation in a sphere of radius 8 / 1 “^ Mpc, as, as 
inferred for each model from the fluctuations in the microwave background detected by COBE (see 
Bunn, Scott, & White 1994 for references). We use the values of as of Bunn et al. (1994) for sCDM 
(cTs = 1.34) and MDM (crs = 0.97), that of Kamionkowski et al. (1994) for CDMo (crs = 0.35), 
and that of Kofman, Gnedin, & Bahcall (1993) for CDMA (crs = 0.79). It should be pointed out 
that Bunn et al. (1994) and Kamionkowski et al. (1994) use the full two-year COBE data, whereas 
Kofman et al. (1993) use the first-year data only. In addition, the fitting procedures of the latter 
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two papers are not as accurate as those of the first.^ However, these errors are comparable to the 
measurement errors and the effect of cosmic variance. 

From £78 and the linear power spectrum, softened at small scales by pressure effects (eq. 4), we 
infer the value of cr/(o = 1) = crjo, the linearly extrapolated fractional gas mass fluctuation in a 
sphere of vanishing radius at the present time. We use the fit to the power spectrum of Bardeen et 
al. (1986) for sCDM, CDMo, and CDMA, and that of Klypin et al. (1993, note added in proof) for 
the cold component in an MDM model. The growth factor is then normalized to D(o = 1) = ajo 
in the ZA (as discussed in §2.2, this also fixes the normalization in the MZA). Figure 1 shows 
D(o) for all four models, assuming a temperature T = 2000(1 + z)K. The models predict very 
different fluctuations on these small scales. Compared to these differences, the uncertainties due to 
the uncertain temperature are fairly small. For our four models, D(o) decreases by 10 — 17% if T is 
increased by a factor of 3. 

The density PDF in the ZA depends only on D, regardless of expansion parameter and 
cosmological model, and therefore the same PDF occurs at different times in different models. This 
is not exactly true for the PDF of r in the ZA, which also depends (weakly) on din D/d In o, because 
of the mapping from (Eulerian) physical space to velocity space. It is also not true for either PDF in 
the MZA, where the growth factor is different from point to point, and does not appear to be simply 
related to the global growth factor of the ZA. However, we have found that all these differences are 
quite small compared to the difference between models, and to the changes with redshift. Thus, the 
PDFs are fairly well characterized by a single number, the growth factor D of the ZA. 

Figure 2 illustrates the change in the density PDF with varying D. For each of the values, 
D =1/8, 1/4, 1/2, and 1 (corresponding to the rms density fluctuation expected in linear theory), 
we show the full density PDF in the ZA (dashed curves), the density PDF in the ZA considering 
only fluid elements expanding along all three principal axes, i.e., for which all three are negative 
(dotted curves), and the PDF in the MZA considering only fluid elements expanding along all three 
principal axes (solid curves). From the latter two sets of curves, it can be seen that, for progressively 
larger values of D, the fluid elements that expand along all three axes occupy a progressively larger 
fraction of the total volume. It is also clear from comparing the dashed and dotted curves that at 
low densities only the fluid elements that expand along all three principal axes are contributing. 

The central result of this paper are the PDFs for p and r at different redshifts, shown in Figs. 
3, 4, and 5, where the different types of curves have the same meaning as in Fig. 2. These figures 
illustrate that the typical density and optical depth in the interstellar medium are substantially 
lower than they would be if the gas were uniformly distributed. In addition, there are dramatic 
differences among models. Particularly low values are seen for the standard CDM model. 

A comparison of the dashed and dotted curves in Figs. 3-5 also shows that the low-density (or 
low-optical-depth) tail of the PDFs comes nearly exclusively from fluid elements which expand along 
all three principal axes. These regions, although initially filling only ~ 8% of the volume, expand so 
much that they occupy most of it by D ~ 2. The similarity of the dashed and dotted curves almost 

^ After the present work had been submitted, Gorski et al. (1994) reported a more accurate fit 
for open models. They obtain = 0.45 — 0.54 for the particular model considered here, somewhat 
higher than the value we used. This makes the results for CDMo even more similar to those for 
CDMA; see below. 
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up to the “knee” of the latter also suggests that it is acceptable to use only the fluid elements with 
Ai < 0 in the determination of the low density parts of the PDFs in the MZA. Of course, this does 
not guarantee the accuracy of the MZA itself, which should be tested by comparison with numerical 
simulations. 


4. DISCUSSION 

In the previous section, we have presented a calculation of the distribution function of the Gunn- 
Peterson optical depth, arising from underdense regions in a photoionized intergalactic medium 
where the gas is moving towards gravitationally collapsed objects. 

The distribution of the optical depth can be directly determined from the observations of 
quasar spectra, in the region where the Lya forest is observed. Some fraction of the spectrum 
is occupied by detectable absorption lines. These should correspond to overdense regions where the 
gas has collapsed and has been shocked. According to the simulations of Cen et al. (1994), the 
properties of the observed lines seem to be consistent with those arising from structures formed 
during gravitational collapse (sheets, filaments and halos) of the photoionized gas. Obviously, our 
analytical treatment cannot give us the distribution of optical depths from regions in the observed 
spectra occupied by such absorption lines. 

In regions outside these absorption lines, one should observe the optical depth corresponding 
to the underdense gas that is left between collapsed structures, and has not been shocked. The 
distribution of optical depths in these regions, once the tails in the Voigt profiles of the strong 
lines have been subtracted, should correspond to the distribution we have calculated in this paper. 
However, we must bear in mind the limitations of our calculation when it is applied to the 
observations. First of all, our calculation of the optical depth caused by a given point of the 
intergalactic medium in Lagrangian space is only an analytical approximation. It is not clear how 
well our approximation of a constant ratio of the eigenvalues of the deformation tensor in underdense 
regions will work in practice. Second, only a fraction /(t) of all the regions yielding a certain optical 
depth T will not have collided with shocked gas, merging into systems yielding the absorption lines, 
and will also not be superposed with absorption from other regions having the same velocity along 
the line of sight. Here, we simply assume this fraction to be unity for low r, but in general /(t) will 
be less than 1, and the fraction of the spectrum having optical depth lower than r will be smaller than 
what we derive. Third, the hydrodynamic effects of the pressure have not been included, except 
for the smoothing of the power spectrum at the Jeans mass in the linear regime, which we have 
calculated under the assumption that T oc . The streaming of the gas away from the centers of 
voids should be slowed down by pressure gradients in the non-shocked gas; the importance of this 
effect was illustrated in the previous section by considering two different temperatures for the gas. 
Finally, thermal broadening will cause some smoothing of the Gunn-Peterson optical depth, and this 
will reduce the fluctuations. 

All of these effects imply that any results obtained from the observed intensity distributions 
using our analysis can only be taken as preliminary indications of the implications that may be 
derived from observations of the Gunn-Peterson effect. Our analytical results should be useful to 
better understand and interpret the numerical results that can be obtained with simulations like in 
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Cen et al. (1994). The simulations are also subject to effects of finite resolution which may cause 
the underdensity of voids to be underestimated when the density gradients become too large. These 
simulations include all of the effects we have mentioned above, and they also predict the number of 
absorption lines at different column densities produced by collapsed regions. In order to match the 
number of predicted absorption lines in a specific theory with the observed ones, the quantity in 
equation (1) has to be fixed, and the amplitude of Gunn-Peterson fluctuations must then agree with 
observations. 

Photoionized gas will first collapse in objects at the Jeans scale. The gas pressure prevents 
the collapse of smaller objects, causing a cutoff on the power spectrum that will cause a prominent 
structure of pancakes surrounding voids (Zel’dovich 1970). The separation between these pancakes 
should be of order the Jeans scale, Xj H[z) = TSTi^^kms"^, where T 4 = T/(10'*K) is the 
temperature of the intergalactic medium, H[z) is the Hubble constant at the epoch when the 
absorption lines are observed, and we have assumed the density in cold dark matter is equal to the 
critical value (e.g., Peebles 1980, eq. [16.5]). As larger scales continue to collapse, the voids become 
more empty and the sheets formed initially turn thinner as they accrete into more massive structures. 
However, because the gas distribution is smooth on the Jeans scale, there is always an intergalactic 
medium pervading all space between collapsed objects. The average separation between absorption 
lines corresponding to collapsed regions can never be smaller than Xj, but can only increase if the 
sheets in voids become so thin that eventually their absorption lines are indistinguishable from the 
general Gunn-Peterson absorption. 

The point we wish to emphasize here is that, as a consequence of the above paragraph, only the 
absorption lines that are, on average, separated by more than Xj can be associated with physically 
distinct structures, if the Lya forest originates indeed from gravitational collapse. Whenever the 
abundance of the lowest column density lines rises above this value, these lines must correspond 
to the fluctuating Gunn-Peterson effect arising from a true intergalactic medium, rather than to 
individual clouds. If they were individual clouds, there is simply not enough mass for them to be 
above the Jeans mass given their abundance, and therefore they could only be pressure confined by 
an intergalactic medium much hotter than the gas in the clouds. One should notice that it is always 
possible to describe a continuum of Lya absorption by a superposition of weak lines, as long as the 
column density distribution and the correlation function of such lines are left as free functions to 
be determined observationally. The fact that such fitting of superposed lines can be done does not 
prove that the lines are associated with real clouds. 

It is interesting to notice that the Gunn-Peterson effect is the only method available to us to 
directly observe underdense matter in the universe. The usual detection of voids from the distribution 
of galaxies (e.g., Kirschner et al. 1981, 1987; Huchra et al. 1983) is done by noticing the absence 
of galaxies on a certain region, and these galaxies are overdense regions on a smaller scale than the 
voids being detected. In general, the primordial fluctuations on scales which have already entered 
the non-linear regime of gravitational collapse cannot be recovered from the distribution of matter 
in the overdense regions, because that is mostly sensitive to non-linear processes of gravitational 
collapse, and the memory of initial conditions is erased. However, the density of matter left in 
underdense regions depends on the initial fluctuations on small scales (e.g., Bernardeau 1992; see 
his Figure 10 for the evolution of the density probability distribution function). If the amplitude 
of fluctuations increases rapidly toward small scales, then the first objects should have collapsed 
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at an earlier epoch, and there has been more time for the matter to flow out of voids. But if the 
power spectrum flattens toward small scales, so that the first objects collapsed more recently, then 
the median density in voids should be higher. Therefore, the determination of the distribution of 
the Gunn-Peterson optical depth, r, for low values of r, can be used to constrain the amplitude 
of primordial fluctuations on small scales, down to the Jeans scale at the temperature expected for 
photoionized gas in the intergalactic medium (T ~ 10'* K). 

The Gunn-Peterson absorption can be observed for two species: HI and Hell. The He I effect, 
occurring at 584 A, is expected to be very small and would be confused with hydrogen lines of lower 
redshift clouds; only a few lines corresponding to high column density systems have been observed 
(Reimers & Vogel 1993). We shall discuss the observations for the other two species. 

4 . 1 . The HI Gunn-Peterson effect 

Most attempts to measure the Gunn-Peterson effect have been based on detecting the discrete 
Lya lines, and then testing if the intensity observed between these lines is consistent with the 
extrapolated continuum of the quasar on the red side of the Lya emission line (e.g., Steidel & Sargent 
1987, Giallongo et al. 1994). Jenkins & Ostriker (1991) presented the distribution of intensities in 
several quasar spectra. They showed that in order to account for the total depression of the flux in 
the region of the Lya forest, there has to be either a uniform Gunn-Peterson effect at the level of 
tgp ~ 0.1 at 2 = 3, or an equivalent absorption due to lines that are too weak and crowded to be 
detected. They also illustrated the effect caused by the presence of a uniform optical depth, r: there 
is a maximum value of the intensity that any pixel in the spectrum can reach, equal to e~'^ relative 
to the extrapolated continuum, except for the instrumental noise. If the absorption between the 
detected lines is to be explained instead by the superposition of weaker, undetected lines, then these 
lines cause not only an additional average depression, but also more fluctuations in the intensity, 
and some pixels are left with very little absorption (see Fig. 3 in Jenkins & Ostriker). 

The higher resolution observations of Webb et al. (1992) of QSO 0000-263, at z = 4.11, clearly 
showed that no uniform Gunn-Peterson optical depth, at the level needed in the models of Jenkins & 
Ostriker where such a continuum absorption was postulated, can be present in the data. Webb et al. 
detected all the absorption lines with Njji > cm“^, and looked at the intensity distribution in 

the regions outside these strong lines. They then modeled this distribution of intensities as arising 
from weaker lines, plus a possible uniform continuum. They found that a steep column density 
distribution, with (3 = 1.7 (where /{Nhi) dNffi oc N^j dNui) cannot be extrapolated to column 
densities Nhi < 10*^cm“^, because that would produce fluctuations of the intensity larger than 
observed. However, if the shallower slope of /3 = 1.3 is assumed for lines with Nhi < 10*3-'!'5cm-2, 
and is extrapolated to zero column density, then both the average depression and the intensity 
fluctuations can be explained, without the need of any uniform absorption. If a continuum absorption 
is present, then the number of weak lines must be even lower, since they must produce a smaller 
average depression, but there must be more of them at Nhi ^ 10*^cm“^ in order to obtain the 
observed intensity fluctuations. 

This turnover in the distribution of column densities in the Lya forest, required by the data of 
Webb et al. (1992), should be expected in a model of gravitational collapse, arising from the lower 
limit to the size of gravitationally collapsed objects, given by the Jeans length (this turnover is seen 
in Figure 3 of Gen et al. 1994 at Nhi — 10*^'®cm“^, although it is probably too pronounced in this 
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case because of the method used to select the lines). In fact, absorption lines with Njji > 10^^cm“^ 
at z = 4 have an average abundance of about one line every 100kms“^ in the observed Lya forest, 
as is easily found from the column density distribution of the observed lines, and the evolution with 
redshift (see, e.g., Petitjean et al. 1993; Murdoch et al. 1986; Lu, Wolfe, & Turnshek 1991; Jenkins 
& Ostriker 1991; Bechtold 1994). This is of the order of the expected Jeans length and, as discussed 
above, these lines must roughly correspond to the weakest ones that can arise from structures that 
have gravitationally collapsed in at least one dimension. Lines with even lower column density should 
form a continuum of absorption, and this true continuum would not be separated into individual lines 
if the resolution and the signal-to-noise were perfect. 

We then conclude that if intergalactic space at 2 ~ 4 is not filled by gas at T ^ 10'* K, heated 
by processes other than photoionization and confining the Lya clouds, then a significant part of the 
fluctuations and the median decrement of the intensity distribution in Figure 2 of Webb et al. (1992) 
must arise from the fluctuating Gunn-Peterson effect. Moreover, the upper limit oi tgp < 0.06 
in Webb et al. is not applicable, because their central assumption of the uniformity of the Gunn- 
Peterson absorption is not valid. 

Following our line of argument, we make the following simple assumption: all the absorption 
lines with Nhi < 10*^cm“^ in the models used by Webb et al. to fit their intensity distribution 
correspond to the fluctuating Gunn-Peterson effect. The distribution of optical depths that these 
lines produce can then be directly compared to our calculations presented in § 3. The first model of 
Webb et al. assumes /3 = 1.7, and a sudden cutoff at Nhi = 10*^cm“^. A uniform Gunn-Peterson 
effect oi Top = 0.04 is then needed. In our interpretation of their results, this represents a lower 
limit to the average optical depth, since in other models incorporating a lower uniform continuum 
the number of lines with Nhi > 10*^cm“^ needs to be reduced, and therefore a larger fraction of 
the total decrement in the intensity must come from weaker lines, which we consider as part of the 
Gunn-Peterson effect. 

Thus, in the second model in Webb et al. , if no lower cutoff is assumed for the column density 
distribution, they find that a value of /3 = 1.3 for Nhi < 10*^'^® cm ^ fits their distribution. The 
average optical depth produced by all these lines is about 0.25 (see Figure 2 of Webb et al. ), and 
approximately 30% of this optical depth is due to lines with Nhi < 10*^cm“^. We then derive the 
higher value tgp = 0.08. 

The turnover of the column density distribution should probably occur smoothly, over a 
wide range in column density, and the optical depth in the Gunn-Peterson effect should then be 
intermediate between the two models discussed above. We therefore adopt a value oi tgp = 0.06 
at z = 4, and assume that this has to be equal to the median optical depth calculated in § 3 for 
different models. 

We start considering our open GDM model, and the GDM model with a cosmological constant. 
Both these models predict the linearly extrapolated rms density fluctuation on the Jeans scale, aj, 
to be about 1 to 1.5 at z = 4. From the solid lines in Figures 4 and 5, we see that the Gunn-Peterson 
optical depth is r < O.lr^ over 25% of the spectrum. If we extrapolate the solid line to higher r, 
assuming it would have a similar shape as the dashed line for the Zel’dovich approximation if regions 
with positive eigenvalues of the deformation tensor were included, 50% of the spectrum should have 
T < 0.2 tu. This should correspond to the observed median Gunn-Peterson optical depth, so we 
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derive ~ 0.3 at z = 4 for these models. On the other hand, for the standard CDM model, the 
fluctuations on the Jeans scale are ~ 3, and the median value of the Gunn-Peterson optical depth 
should be only r ~ O.OItu, implying the much larger value ~ 6 . For models with lower amplitude 
of fluctuations at the Jeans scale, like the MDM model, the required value of is lower. 

In terms of the baryon density of the universe fij, and the intensity of the ionizing background, 
Jhi ,-21 (defined as in eq. [4] of Miralda-Escude & Ostriker 1992), the optical depth in equation 
( 1 ) would be given as: 


Qih ? \^ ^ / i + ^ \® 

0.015y /i \ 5 y H{z) Jhi ,-21 

where T 4 is the gas temperature in units of 10'* K, H[z) is the Hubble constant at redshift z, and we 
have assumed a helium abundance Y = 0.24, and that the helium is doubly ionized. For the open 
CDM model we use here, H[z)/Ho = 7.4 at z = 4, and for the model with a cosmological constant, 
H[z)/Ho = 7.1. Fixing fij to the value required by primordial nucleosynthesis, we need Jhi ,-21 — 1 
for these two models, in good agreement with determinations of this intensity from the proximity 
effect (Bajtlik, Duncan & Ostriker 1988; Lu, Wolfe, & Turnshek 1991; Bechtold 1994). On the other 
hand, for the standard CDM model, we need Jhi ,-21 — 0 . 2 . 

The simulations of Cen et al. (1994) were done for the same CDM model with cosmological 
constant that we have used here. According to their results, the number of absorption lines is 
correctly predicted if Jhi ,-21 — 1; which is the same value that we find is needed to explain the 
Gunn-Peterson effect. Thus, this model seems to be satisfactory, and probably any models with rms 
density fluctuations on the Jeans scale near unity should give similar results. 

If the density fluctuations are larger, like in standard CDM with the COBE normalization, then 
systems on larger scales will have collapsed, which will have higher total column densities of gas. 
The column densities of neutral hydrogen will be even more increased relative to models with smaller 
density fluctuations, because a lower intensity of the ionizing background is required to explain the 
amplitude of the Gunn-Peterson effect. Thus, the turnover in the column density distribution would 
occur at a column density that is too high. Such models would therefore be in conflict with the 
observations. Their basic difficulty is that there is too much gas in collapsed systems on large scales, 
and too little left in the intergalactic medium, so that the contrast between the high column density 
lines and the Gunn-Peterson effect is too large. 

At the opposite extreme, if the density fluctuations are too small, there should be too few 
collapsed systems producing high column density absorption lines, and the Gunn-Peterson effect 
would be too large, or equivalently, there would be too many superposed weak lines, giving a 
very steep column density distribution. Models of this type, such as the MDM model, also have 
difficulty in forming a sufficiently large number of halos where the gas can dissipate and form 
damped Lya systems and high redshift quasars with the observed numbers (Mo & Miralda-Escude 
1994; Kauffmann & Gharlot 1994; Ma & Bertschinger 1994). 

We then conclude that models with an rms density fluctuation of order unity at the Jeans scale 
are favored by the observations of the Lya forest, and by the distribution of intensities in Webb et 
al. (1992). As we have discussed before, this is only a preliminary result, which needs to be further 
investigated with numerical simulations. However, it seems clear that the ratio of the Gunn-Peterson 
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optical depth to the number of individual absorption lines will constrain the primordial rms density 
fluctuations as we have discussed in this paper. 


4 . 2 . The Hell Gunn-Peterson effect 


For the He II Gunn-Peterson effect, all our equations are equally applicable as long as most of the 
helium in the intergalactic medium is twice ionized. In this case, the density of singly ionized helium 
is also approximately proportional to as discussed in § 2, since the recombination coefficient 

for Hell depends on temperature almost in the same way as the HI recombination coefficient. The 
abundance of He I is always very small, and has a negligible effect on the equilibrium of He II. 
Equation (1) is changed to 


Tu,HeII — 2.5 X 10^ 


0.5 

0.015 IT 


I + 2V Ho 


(27) 


Jakobsen et al. (1994) have recently reported the first observation of the spectrum of a quasar 
at a rest-frame wavelength immediately below 304 A. A cutoff in the flux was observed at this 
wavelength, and interpreted as the Hell Gunn-Peterson effect. The trough is observed at 2 ~ 3, 
and the flux decrement is more than a factor of 5 at the two-sigma level. 

From our discussion above, it is clear that in a theory of gravitational collapse for the Lya 
forest, such a trough must be due to the Gunn-Peterson effect, since any weak lines causing it would 
have to be overlapping. However, because of the fluctuations in the Gunn-Peterson effect, we expect 
that some regions in the spectrum should have low optical depth, and show up as “gaps” in the 
He II trough, which could be observed in a spectrum of better resolution and signal-to-noise than in 
Jakobsen et al. (1994). Such gaps should be a very sensitive probe to the most underdense voids in 
the intergalactic medium. 

In general, the lower limit required on the fraction of singly ionized helium in the intergalactic 
medium to give a certain He IIGunn-Peterson optical depth is larger than in the case of a uniform 
medium, because the density of the gas is much lower than average in the voids. However, the 
constraints on the spectrum of the ionizing background are not modified, since this depends only on 
the ratio of the H land He IIGunn-Peterson effects. 


5. CONCLUSIONS 


In this paper, we have considered the model of gravitational collapse from primordial density 
fluctuations for the Lya forest, where photoionized gas in the intergalactic medium collapses into 
the structure forming on scales above the Jeans mass. The observations of the Gunn-Peterson 
effect and the low column density lines in the Lya forest are unique among all other cosmological 
observations, in the sense that they allow us to directly observe matter in underdense regions, as 
well as moderately overdense matter in the infall regions around collapsed objects. An analytical 
approximation was used to calculate the distribution of the density field and the Gunn-Peterson 
optical depth from regions that are initially expanding in all directions in comoving coordinates. 
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We have shown that, in this theory, the distribution of intensities observed by Webb et al. 
(1992) in the 2 = 4.1 quasar QSO 0000-263 gives evidence that the fluctuating HI Gunn-Peterson 
effect has been detected. Over most of the spectrum of a quasar outside strong lines the optical 
depth originates from underdense regions, and its typical value, which is found to he tgp — 0.06 at 
2 = 4, is therefore significantly lower than the value expected for a uniform intergalactic medium 
with the baryon density implied by primordial nucleosynthesis. 

Models for structure formation in which the fractional rms density fluctuation is close to unity 
at the Jeans scale at z = 4 imply a value of that is consistent with the measured intensity 
of the ionizing background from the proximity effect, Jhi — 10“^^ erg s“^ cm“^ Hz~^ sr“^. From 
the results of Cen et al. (1994), these models are also consistent with the observed column density 
distribution of the absorption lines, for the same value of Jhi- A much higher or much lower 
value for the density fluctuations should lead to different predictions for the relative strength of 
the fluctuating Gunn-Peterson effect and the individual absorption lines. These conclusions are at 
this point preliminary, but further studies with numerical simulations for several models, and better 
observations of high redshift quasar spectra at high resolution, should give more definite conclusions. 

There are three complications that must be considered, which can affect the properties of 
the intergalactic medium in ways that cannot be incorporated in numerical simulations using only 
basic physical principles. The first is the fact that the intensity and the spectrum of the ionizing 
background will fluctuate due to random variations in the number of emitting sources, as well as 
the number of absorbers (Zuo 1992a,b; Fardal & Shull 1993). The second is that the temperature 
of the intergalactic medium should also fluctuate, and may not be a function of the density only, 
because it depends on the initial temperature to which the gas is heated at reionization. This initial 
temperature can vary substantially from place to place, depending on the spectrum of the sources 
reionizing the medium (Miralda-Escude & Rees 1994). In addition, the spectrum of the ionizing 
background should also vary depending on the random fluctuations in the number of sources and 
absorbers, and this affects the equilibrium temperature of the gas. Finally, energetic events in the 
collapsed objects, such as supernova explosions and quasars, can eject gas into the intergalactic 
medium and heat it to high temperatures. 

The latter process is the one that can change more drastically the properties of the Gunn- 
Peterson effect. If ejection of gas is important, then the voids can be refilled and more absorption 
lines can be produced by clouds not associated with gravitational collapse. Such lines could constitute 
a different population of objects, and should have smaller transverse sizes than the gravitationally 
collapsing systems. The Gunn-Peterson effect could decrease much more slowly with time if gas were 
constantly ejected into voids. Recent observations of GIV in systems with Nhi ~ 3 x 10^'*cm“^ 
(Gowie et al. 1995; Fan & Tytler 1995) suggest that the metallicity is not very low compared to 
systems of higher column density. If this is also true for much lower column density lines in the 
Lya forest, then the intergalactic medium must have been contaminated with metal-rich, ejected gas 
at some point in the past. This could have occurred at a very early epoch, when the intergalactic 
medium was being reionized. In that case, the metals could have been ejected from systems formed 
from neutral gas on very small scales, and the ejection velocities might be too small to significantly 
affect the evolution of the voids in the intergalactic medium. However, if the metal enrichment 
was from systems that can eject gas at velocities comparable to those generated by the collapsing 
structures, then the assumptions we have made here would be invalidated. This could then allow 
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models with high initial density fluctuations to be made consistent with observations of the Gunn- 
Peterson effect. 

If it is found that the simple case where all these complicated effects do not have a strong 
influence on the intergalactic gas can predict the observed properties of the Lya forest, then it 
should be clear that most of the volume in the intergalactic medium at z = 4 should not have been 
affected by shock waves from explosions. Fluctuations in the gas temperature and the intensity of the 
background cannot drastically alter the median Gunn-Peterson optical depth. The amplitude of the 
Gunn-Peterson fluctuations should then provide a strong constraint on the primordial fluctuations 
at small scales. 
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FIGURE CAPTIONS 


Fig. 1 : Fractional rms density fluctuations at the Jeans scale, aj, calculated using linear 
perturbation theory, as a function of the expansion factor, a = (1 + 2 :)“^, in the four models of 
structure formation considered in this paper. The assumed baryon temperature is 2000(1 + z)K. [aj 
also corresponds to the growth factor D[a) in the Zel’dovich approximation). 

Fig. 2 : Cumulative probability distribution function (PDF) of the density (i.e., fraction of the 
volume at densities lower than a certain value p) as a function of \og[p/p), where p is the average 
density, for different values of the fractional rms density fluctuation calculated in linear theory, 
a = 1/8, 1/4, 1/2, and 1. The dashed lines were obtained in the Zel’dovich approximation (ZA), 
considering all matter elements. The dotted lines also correspond to the ZA, but considering only 
those matter elements which expand along all three principal axes. The solid lines were obtained 
in our modified Zel’dovich approximation (MZA) using a flat, CDM-dominated cosmological model, 
and again considering only matter elements which expand along all three principal axes. 

Fig. 3 : Cumulative PDFs of the density for four models of structure formation described in the text. 
The meaning of the line types is the same as in Fig. 2. Within each type, the four lines represent 
(from right to left) redshifts z = 4,2, 1, and 0. The assumed gas temperature is T = 2000(1 + z) K. 

Fig. 4 : Cumulative PDFs of the optical depth (i.e., fraction of a quasar spectrum sampling matter 
of an optical depth lower than a given value r), as a function of log(T/Tu), where is the optical 
depth expected if the matter were uniformly distributed and expanding with the Hubble flow. The 
models and the redshifts are the same as in Fig. 3 (but note that the lower redshift lines are off 
the scale in some of the plots). Line types are as in Fig. 2. The assumed gas temperature is 
T = 2000(1 + 2) K. 

Fig. 5: The same as Fig. 4, except at a gas temperature T = 6000(1 + z)K, which increases the 
Jeans scale by 3^/^ sa 1.7, and therefore lowers the amplitude of the fluctuations. 



